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Maximally localized Wannier functions (MLWFs) based on Kohn-Sham band-structures provide a 
systematic way to construct realistic, materials specific tight-binding models for further theoretical 
analysis. Here, we construct MLWFs for the Mn Cg bands in LaMnOa, and we monitor changes in 
the MLWF matrix elements induced by different magnetic configurations and structural distortions. 
From this we obtain values for the local Jahn- Teller and Hund's rule coupling strength, the hopping 
amplitudes between all nearest and further neighbors, and the corresponding reduction due to the 
GdFeOa-type distortion. By comparing our results with commonly used model Hamiltonians for 
manganites, where electrons can hop between two "eg-like" orbitals located on each Mn site, we 
find that the most crucial limitation of such models stems from neglecting changes in the underlying 
Mn(ci)-0(p) hybridization. 



I. INTRODUCTION 

The theoretical description of complex transition metal 
oxides and similar materials is very often based on effec- 
tive tight-binding (TB) models, i.e. a representation of 
the electronic structure within a certain energy region in 
terms of localized atomic-like orbitals. Simple TB mod- 
els with a small number of orbitals can be used to study 
the essential mechanisms governing complex physical be- 
havior, such as for example that found in the colossal 
magneto-resistive manganites. 

The electronic properties of manganites Ri-xAxM-iiO^ 
{R: trivalent rare earth cation, A: divalent alkaline 
earth cation) are often described within an effective "two- 
band" TB model, where electrons can hop between the 
two eg levels on each Mn site. The corresponding Hamil- 
tonian typically also contains several local terms describ- 
ing the coupling of the Cg states to the t2g core spin, 
to the Jahn- Teller (JT) distortion of the oxygen octa- 
hedra, and/or the electron-electron Coulomb repulsion. 
It has recently been shown, that such a model (with pa- 
rameters obtained partly from first principles calculations 
and partly by fitting to experimental data) is able to re- 
produce the basic structure of the phase diagram as a 
function of doping and temperature found in manganite 
systems such as Lai_a;(Ca,Sr)a:Mn03.^ 

An elegant and systematic way to obtain realistic 
(materials-specific) TB models is the construction of 
maximally localized Wannier functions (MLFWs) from 
the Kohn-Sham states calculated using density functional 
theory (DFT).Sl DFT calculations are known to give a 
realistic description of electronic structure for systems 
where electronic correlation effects are not too strong 
Furthermore, for materials where correlation effects are 
important, a Wannier representation of the Kohn-Sham 
band structure can be used to define a subset of orbitals 
(the "correlated subspace"), which can then be used as 
basis for a more elaborate treatment of correlation ef- 
fects beyond standard DFT. This is done for example in 
DFT+DMFT (DMFT = dynamical mean-field theory) 



calculations,™^ which aim at an accurate quantitative 
description of materials where electronic correlation can- 
not be ignored. 

In this work we construct MLWFs corresponding to the 
Mn Cg states for LaMnOa, the parent compound for many 
manganite systems, based on DFT calculations within 
the generalized gradient approximation (GGA). We cal- 
culate the real space Hamiltonian matrix elements in the 
MLWF basis for different structural modifications and 
for different magnetic configurations, and we compare 
the obtained results with assumptions made in commonly 
used two band TB models. 

Ovir analysis is closely related to earlier work presented 
in Ref. ,10i which examined the validity of the two band 
picture by fitting TB model parameters (including the 
hopping between nearest and next-nearest neighbors) to 
the DFT band structure obtained within the local density 
approximation (LDA). The approach based on MLWFs 
used in the present work is less biased and more generally 
applicable, and thus allows for a more systematic anal- 
ysis than the manual fitting of TB parameters discussed 
in Ref. [TOl It is also well suited for the construction 
of the correlated orbital subspace used for DFT+DMFT 
calculations .'2' 

This paper is organized as follows. In the following 
section we describe the theoretical background of our 
work. Thereby, Sec. II A summarizes the effective two 



band model that is often used for a theoretical treatment 
of manganites, Sec. |IIB| presents the definition of the ML- 
WFs, Sec. |II C| describes the various structural modifica- 
tions of LaMnOa investigated throughout this work, and 
Sec. |IID| lists some of the calculational details. The pre- 
sentation of results starts with the case of the ideal cubic 
perovskite structure in Sec. |III A[ The individual effects 
of the staggered JT and the GdFeOs-type distortions are 
then presented in Sees . |III B] and [iTl C | respectively. This 
is followed by the results for the combined distortion in 
Sec. IIIIDI and the construction of a refined TB model 



and its application to the full experimental structure of 
LaMnOs in Sec. |IIIEl Finally, the most important results 
and conclusions are summarized in Sec. II VI 
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FIG. 1: (Color online) Different structural modifications of 
LaMnOa investigated in this work, viewed along the [001] 
direction: (i) ideal cubic perovskite, (ii) purely Jahn- Teller 
distorted, (iii) purely GdFeOa-type distorted, and (v) exper- 
imental Pbnm structure. Pictures have been generated using 
VESTA.IIB 



II. METHOD AND THEORETICAL 
BACKGROUND 

A. Effective two-band models for LaMnOg 

LaMnOa crystallizes in an orthorhombically dis- 
torted perovskite structure with Pbnm space group (see 
Fig.[T|ir), and A- type antiferromagnetic (A -AFM ) order of 
the magnetic moments of the Mn cations-^i^ES! The devi- 
ation from the simple cubic perovskite structure (Fig. [TJ) 
can be decomposed into a staggered JT distortion of 
the MnOg octahedra within the x-y plane (Fig. [Tji), the 
so-called GdFeOs-type (GFO) distortion, consisting of 
collective tiltings and rotations of the oxygen octahedra 
(Fig. [Tjii), and "the rest", i.e. displacements of the La 
cations from their ideal positions plus a homogeneous or- 
thorhombic strain (Fig. 1}^).^° 

The electronic structure of LaMnOa close to the Fermi 
energy is dominated by Mn 3d states, which are split by 
the cubic component of the crystal field into the lower- 
lying three-fold degenerat e t2a an d the higher-lying two- 
fold degenerate eg states.'^^^^^^ The formal electronic 
configuration Mn'^"'' : [Ar] 3d* leads to a high spin state of 
the Mn cation with fully occupied local majority spin t2g 
states and one electron per local majority spin eg state, 
while both t2g and eg minority spin states are empty. 

Based on this electronic structure, the theoretical de- 
scription of manganites often employs an effective two- 
band TB picture, where electrons can hop between the 
two eg levels on each Mn site. This hopping is facilitated 
by hybridization with the oxygen 2p states, which, how- 
ever, are not explicitly included in the TB model. It is 



therefore understood, that the "atomic" eg states used 
in the TB model are indeed somewhat extended Wan- 
nier orbitals that also include the hybridization with the 
O 2p states. In contrast, the three t2g electrons are as- 
sumed to be tightly bound to a specific Mn site where 
they give rise to a local "core spin" S = 3/2. This core 
spin then interacts with the valence eg electron spin via 
Hund's rule coupling. In addition, a JT distortion of the 
surrounding oxygen octahedron splits the two eg levels 
on the corresponding Mn site, whereas elastic coupling 
between neighboring oxygen octahedra gives rise to a co- 
operative effect. The GFO distortion in this picture is 
usually assumed to simply reduce the effective hopping 
amplitudes between neighboring Mn sites due to the re- 
sulting non-ideal Mn-O-Mn bond angle. In addition, a 
local electron-electron interaction between electrons oc- 
cupying the same Mn site can be included in the model.'^ 

The electronic Hamiltonian for such a model can be 
expressed as: 



H — i/kin + H] 



local 



(1) 



where 



4(R-|-AR)o-'^aRcr + h.C. (2) 

a,b,'R,AIl,a 

describes the electron hopping between orbital \a) (spin 
a) at site R and orbital \b) at site R -I- AR, and it is 
assumed that all sites are translationally equivalent, so 
that the hopping amplitudes iab(AR) depend only on the 
relative position between the two sites. 

Representing the eg orbital subspace within the usual 
basis |1) = |3z^ — r^) and |2) — \x'^ — y^), and assuming 
cubic symmetry, the nearest neighbor hopping along the 
three cartesian directions has the following form: 



t(±acz) 



t(±acx) = t 



t(±acy) 




(3) 
(4) 

(5) 



Here, is the lattice constant of the underlying cubic 
perovskite structure. The hopping t' between two neigh- 
boring \x'^ — ?/^)-type orbitals along z is small due to 
the planar shape of this orbital, and it is therefore often 
neglected. In this case, the nearest neighbor hopping de- 
pends only on a single parameter t, the hopping along z 
between |3z^ — r^)-type orbitals. 

H\oced contains all local interaction terms included in 
the model, i.e. Hund's rule coupling with the t2g core 
spin, the JT coupling to the oxygen octahedra distortion, 
and eventually also the electron-electron interaction. In 
this work we will discuss only the Hund's rule and JT 
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coupling, which are of the form: 

Hnund = -J^^R - Sr, , and 



(6) 



R 



H 



JT 



R,.(T,a,b 



(7) 

Here, J is the Hund's rule coupling strength and Sr, is 
the t2g core spin at site R, which in the following we 
will consider as classical vector normalized to |Sr| = 1. 

= J2a,a.cr' clu^Tc^a' CaRa' IS the Corresponding eg va- 
lence spin, A describes the strength of the JT coupling, 
and To-ct' are the usual Pauli matrices. The quantities 
and describe the JT distortion of the oxygen 
octahedron surrounding site R: 



Q 



1 



R 



2^/2 



(dR-dl) 



(8) 



(9) 



where df^, d^, and d^ indicate the 0-0 distances along 
the X, y, and z directions, corresponding to the oxygen 
octahedron located at site R. 



B. Mciximally localized Wannier functions 

As is well known from basic solid state physics, the 
eigenfunctions within a periodic crystal potential are ex- 
tended Bloch waves, classified by a wave- vector k and 
a band-index m. These Bloch waves can in turn be ex- 
pressed as a Bloch sum of "atomic-like" localized TB ba- 
sis functions or Wannier functions. For an isolated group 
of N Bloch states IV'mk) , i-e. a group of bands that are en- 
ergetically separated from all lower- or higher-lying bands 
throughout the entire Brillouin zone (BZ), a set of N lo- 
calized Wannier functions |w„t), associated with lattice 
vector T, is defined via the following transformationPE^l 



\WnT) 



V 

(2^ 



BZ 



N 

E 



ik-T 



dk. (10) 



Here, U^'^^ is a unitary matrix mixing Bloch states at 
wave- vector k. Different U''^'' lead to different Wannier 



orbitals, which are not uniquely determined by Eq. ( 10 1 



However, Marzari and Vanderbilt showed that a unique 
set of maximally localized Wannier Junctions (MLWFs) 
can be obtained by minimizing the total quadratic spread 
of the Wannier orbitals, defined asP 



N 



(11) 



where (O), 



{Wn0\O\w„ 



For the case of entangled Bloch bands, i.e. bands 
that are not energetically separated from other groups 
of higher- or lower- lying states, an energy window 

[Erain, ^^max] Can be defined such that there are n'^^^ > N 
Bloch bands within the energy window at each k vec- 
tor, and then an A'^-dimensional manifold of mixed Bloch 
states is obtained 



1^; 



dis \ 



E 



u: 



dis(k) 
lin 



IV-ik) 



(12) 



The corresponding Wannier functions can then be ob- 
tained from the mixe d Bloch states by replacing jf/'mk) 



with |-0;J^k> in Eq. ( 10 ) . The unitary rectangular N]^^l x N 



matrix U'^"^'^''^ is also uniquely determined by the condi- 
tion of maximal localization, i.e. it can be obtained by 
minimizing ^ }^ 

Once a set of MLWFs is obtained, the corresponding 
Hamilton matrix, H*-^^ (k), is constructed by a unitary 
transformation : 

H(w) (k) = (u^'^y (U'l'^f'^yH^^) (k) U'l'^O^^U^'') , 

(13) 

from the (diagonal) Hamilton matrix in the Bloch basis, 

Hnm (k) = e„k^nm, with eigenvalues e„k- The MLWF 
Hamiltonian in real space is then calculated as a Fourier 
transform of H*^^^ (k), which in practice is replaced by a 
sum over Nk points in k-space: 



k 



ik-T 



(14) 



Thus, the real space representation of the Hamiltonian 
in the MLWF basis is equivalent to a TB description 
of the full Hamiltonian within the corresponding orbital 
subspace: 



H = 



T,AT 



nT+AxCmT 



h.c. 



(15) 



where CmT is the annihilation operator for an electron in 
orbital |u>,„t)- The real space MLWF matrix elements 
^Tim can therefore be interpreted as hopping amplitudes 
within a TB picture of MLWFs [compare Eq. ( 15 1 with 
Eq. Note that AT in Eq. ([l5]) refers to lattice vec- 
tors, whereas AR in Eq. p| refers to Mn sites. The 
subscripts n and m in Eq. ( [T5| can thus in general indi- 
cate both site and orbital/spin character (for cases with 
more than one site per unit cell). 

For the case when MLWFs are constructed from an 
isolated set of bands, the TB model, Eq. (15), exactly 



reproduces the band dispersion within the correspond- 
ing energy window. For the entangled case, the energy 
bands calculated from Eq. (15) do not necessarily have 



to coincide with the underlying Bloch bands. 
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C. Structural decomposition 

To analyze the effect of the various distinct structural 
distortions within the experimental Pbnm structure on 
the electronic properties of LaMnOa we investigate sev- 
eral different atomic configurations (similar to Ref. [rUjl : 

(i) The ideal cubic perovskite structure (Fig. 

(ii) A purely JT distorted structure (Fig. [iji), which 
results from alternating long and short 0-0 dis- 
tances within the x-y plane, i.e. a staggered JT 
distortion Qi^ = ±(5§ and = 0. This distortion 
doubles the unit cell within the x-y plane, leading 
to new in-plane lattice vectors aji = ac(y + x) and 
bii = ac(y — x) and tetragonal symmetry. 

(iii) A purely GFO-distorted structure (Fig. [l}ii), re- 
sulting from rotations of the oxygen octahedra 
around the z direction and octahedral tilts away 
from z, alternating along all three cartesian di- 
rections. This distortion quadruples the unit cell 
compared to the undistorted structure (i), yield- 
ing orthorhombic Pbnm symmetry. The result- 
ing in-plane lattice vectors are identical to those 
of structure (ii) and the new lattice vector along z 
is Ciii = 2acZ. 

(iv) A superposition of JT and GFO distortion, which 
also leads to orthorhombic Pbnm symmetry and 
unit cell vectors unchanged with respect to struc- 
ture (iii). 

(v) The full experimental structure (Fig. 
with orthorhombically strained lattice vectors 
(|av| 7^ |bv| |cv|, resulting in ^ 0) and 
displaced La cations compared to structure (iv) . 

For each of these structural modifications we use the 
same volume V = 60.91 per formula unit as in the 
experimentally observed Pbnm structure.!^ This leads 
to a cubic lattice parameter Oc — 3.9345 A, which devi- 
ates only by 0.8% from the value we obtain by volume 
optimization for the ideal perovskite structure within 
GGA. For the positions of the O anions in structures 
(ii) and (iii) we use the same decomposition of structure 
(iv) into pure JT and GFO components as described in 
Ref. [in] (see Table For the cases with A-AFM order, 
the unit cell is doubled in z direction for both (i) and (ii) 
structures in order to accommodate the magnetic order, 
thus changing the symmetry to tetragonal in case (i). 

Starting from the ideal cubic perovskite structure, we 
analyze the effect of a specific distortion by gradually 
increasing the amount of this distortion, i.e. we perform 
series of calculations using a linear superposition of the 
Wyckoff positions in the cubic perovskite structure and 
in structure (x): 

R(a,) = (l-a,)R« + a,R("), (16) 



TABLE L Wyckoff parameters of the 0(4c), (x, y, 0.25), 
0(8d), (x, y, z), and La(4c), {x, y, 0.25), sites for the vari- 
ous structural configurations used in this work (compare with 
Table I in Ref. fTO)l . 







Expt. (ReL [T7j) 


(ii) 


(iii) 


(iv) 


(v) 


0(4c) 


X 


-0.0733 


0.0 


-0.0733 


-0.0733 


-0.0733 




y 


-0.0107 


0.0 


-0.0107 


-0.0107 


-0.0107 


0(8d) 


X 


0.2257 


0.2636 


0.2121 


0.2257 


0.2257 




y 


0.3014 


0.2636 


0.2879 


0.3014 


0.3014 




z 


0.0385 


0.0 


0.0385 


0.0385 


0.0385 


La(4c) 


X 


0.0063 


0.0 


0.0 


0.0 


0.0063 




y 


0.5436 


0.5 


0.5 


0.5 


0.5435 



and vary ax between and 1. The following cases are 
considered: (x =ii) (pure JT distortion), (x =iii) (pure 
GFO distortion) , {x =iv) (combined JT and GFO dis- 
tortions) . 

D. Computational details 

We perform spin-polarized first principles DFT calcu- 
lations using the Quantum-ESPRESSO program pack- 
g^ggliSJ |.]^g GGA exchange-correlation functional of 
Perdew, Burke, and Ernzerhof,'!^ and Vanderbilt ul- 
trasoft pseudopotential^^Sl in which the La (5s, 5p) and 
Mn (3s, 3p) semicore states are included in the valence. 

Convergence has been tested for the total energy 
and total magnetization using the ideal cubic perovskite 
structure and ferromagnetic (FM) order. We find the to- 
tal energy converged to an accuracy better than 1 mRy 
and the total magnetization converged to an accuracy 
of 0.05 /iB for a a plane- wave energy cut-off of 35 Ry 
and a F-centered 10 x 10 x 10 k-point grid using a Gaus- 
sian broadening of 0.01 Ry. These values for plane- 
wave cutoff and Gaussian broadening are used in all 
calculations presented in this work. The 10x10x10 k- 
point grid is used in all calculations for the cubic struc- 
ture (i), whereas appropriately reduced k-point grids of 
10x10x5, 7x7x10, and 7x7x5 are used for the struc- 
tures with unit cell doubled in the z direction, doubled 
in the x-y plane, and quadrupled, respectively. 

After obtaining the DFT Bloch bands within GGA, 
we construct MLWFs using the wannier90 program inte- 
grated into the Quantum-ESPRESSO packagct^U Start- 
ing from an initial projection of the Bloch bands onto 
atomic d basis functions |3z^ — r^) and \x'^ — y^) centered 
at the different Mn sites within the unit cell, we obtain 
a set of two eg-like MLWFs per spin channel for each 
site. The spread functional (both gauge-invariant and 
non-gauge-invariant parts) is considered to be converged 
if the corresponding fractional change between two suc- 
cessive iterations is smaller than 10~^°. For cases with 
entangled bands a suitable energy window is chosen as 
described in the corresponding "Results" section. 
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FIG. 2: (Color online) Projected DOS and band structure 
along high symmetry lines within the BZ calculated for the cu- 
bic structure (i) and both FM and A-AFM order. Filled (red) 
areas and solid (green) lines represent the projected DOS cor- 
responding to Mn(eg) and Mn(f2g) states, respectively, while 
dashed (blue) lines represent the site and orbitally averaged 
projected DOS corresponding to the O p states. For the A- 
AFM case the left (right) panel corresponds to local majority 
(minority) spin projection. In the band structure plots, the 
dispersion calculated from the eg-like MLWFs are represented 
by thick (red) lines. The Fermi level is indicated by the hori- 
zontal dashed lines. 



III. RESULTS AND DISCUSSION 
A. Perfect cubic perovskite — structure (i) 

The projected densities of states (DOS) and band 
structure calculated for LaMnOa in the ideal cubic per- 
ovskite struct ure ( i) for both FM and A-AFM order are 
shown in Fig. [2p2J 

A metaUic state is obtained for both FM and A- 
AFM order, in agre ement with previous band-structure 
calculations.™™^ The projected DOS show that the 
(local) majority spin bands around the Fermi energy have 
mainly Mn(eg) character and are half-filled while the (lo- 
cal) minority spin bands with mainly Mn(eg) character 
are unoccupied, as expected from the formal electron con- 
figuration. Bands with Mn(t2g) character are lying just 
below the Mn(eg) bands, and slightly overlap with the 




FIG. 3: (Color online) Real space representation of the ML- 
WFs for majority and minority spin projections in the cubic 
structure (i) with FM order, projected on the x-z plane pass- 
ing through Mn (large/blue sphere) and O (small/red spheres) 
sites (in arbitrary units). 



latter. 0(p) bands are located below the Mn(i2g) bands 
(between ~ 6-12 eV) and are fully occupied. The strong 
hybridization between Mn(d) and 0(p) electrons can be 
seen from the substantial amount of Mn((i) character in 
the energy range around 8 eV, i.e. towards the bottom of 
the bands with predominant 0(p) character. The states 
above the Mn(eg) bands have predominantly La((i) char- 
acter. 

One can see from the band structures depicted in Fig. [2] 
that for the FM majority spin channel the bands with 
predominant eg character are nearly completely isolated 
from both higher and lower-lying bands, while for the FM 
minority spin channel and in the A-AFM case, the "cg 
bands" overlap strongly with other bands (with mostly 
Mn(t2g) minority and La(c?) character). As described 
in section IIIDl in order to construct e„-like MLWFs for 



the various cases, we define an energy window for the 
disentanglement procedure [see Eq. (12)], and then ini- 



tialize the Wannier functions from a projection of the 
Kohn-Sham states within that energy window on atomic 
Eg wave- functions (see Ref . [1^ . A suitable energy win- 
dow is chosen based on the eg projected DOS and cal- 
culated band structure (see discussion below for more 
details) . Two MLWFs per spin channel for the single Mn 
site within the cubic unit cell are constructed for FM or- 
der, and two pairs of MLWFs, localized at the two Mn 
sites within the magnetic unit cell, are constructed for 
A-AFM order (for global spin up projection only). 

Figure [3] shows the real space representation of the 
two eg-like MLWFs for both majority and minority 
spin and FM order, calculated for an energy window 
of [12.0, 17.0] eV and [15.9, 20.0] eV, respectively. The 
shape of the MLWFs resembles the antibonding cr* char- 
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acter of hybridization between Mn(eg) and 0{p) states in 
this energy range. The hybridization is notably stronger 
for the majority spin MLWFs (individual spread per WF 
2.90 compared to 1.65 for the minority spin ML- 
WFs), which is due to the smaller energy separation be- 
tween the atomic Mn(eg) and 0{p) levels for the majority 
spin channel. The difference between the real space rep- 
resentation of the MLWFs for FM and A-AFM order (not 
shown here) is more subtle. A quantitative comparison 
of the corresponding differences in the real space Hamil- 
ton matrix elements between MLWFs will be presented 
below. 

The dispersion calculated from the obtained Cg ML- 
WFs is also shown in Fig. |2] It can be seen that even in 
the cases with strongly entangled Cg bands (FM minor- 
ity spin and A-AFM) the MLWF bands follow certain 
DFT bands almost exactly, except around some band 
crossings with higher lying La d bands. This represents 
the fact that within cubic symmetry the eg states can- 
not hybridize with the t2g bands, and hybridize only very 
weakly with the La d states. 

In order to reproduce the two majority spin bands 
around the Fermi energy for the FM case, the lower 
bound of the energy window, i?mini has to be above 
the lower peaks in the Mn(eg) projected DOS at around 
10.5 eV and 8 eV, which correspond to the bonding com- 
bination of hybridized atomic 0{p) and Mn(eg) states. If 
these bands are included in the energy window, the bond- 
ing and antibonding combinations of atomic orbitals be- 
come disentangled and the eg Wannier functions become 
essentially "atomic-like" (compare also with the case of 
SrVOa described in Ref. "^j. On the other hand, varying 
-E'min within 0.4 eV below the F point energy of the Cg- 
like bands changes the MLWF bands by less than 1 meV 
for any k. Similarly, varying the upper bound of the 
energy window has only minor influence on the result- 
ing MLWF bands, due to the negligible hybridization of 
the Cg states with higher-lying bands. Additional test 
calculations for different k-point grids showed that the 
MLWF band structure is converged within 0.5 meV at 
any k-point for the 10 x 10 x 10 grid which was used for 
the energy window test calculations. 

We now turn to the analysis of the hopping parame- 



ters, i.e. the real space matrix elements h^^, Eq. (|14| 



between MLWFs located at different Mn sites. The mag- 
nitudes of all calculated hopping parameters for the FM 
majority spin case are shown in Fig. [4] It is noticeable 
that the hopping amplitudes along the three cartesian 
axes are most dominant and that their decay with dis- 
tance is rather slow, so that the terms corresponding to 
inter-site distances of 2ac and Soc are of comparable mag- 
nitude as the hopping between next-nearest neighbors for 
which |AR| = V2ac. 

The exact MLWF representation in terms of H(^)(k) 
is well suited for further numerical calculations, e.g. 
within a DFT+DMFT approach. On the other hand, 
for the analysis of specific physical mechanisms within 
a semi-analytical TB model, one generally wants to use 




10' 1 2"- 3"= 2 3 4 ?* 



lARI/a, 

FIG. 4: (Color online) Magnitude of all calculated non-zero 
hopping parameters for FM order in the ideal cubic structure 
as a function of the inter-site distance |AR| (open circles: 
hopping along the unit cell directions; open diamonds: hop- 
ping between next-nearest neighbors; filled circles: all other 
hoppings). Inset: Comparison of the full MLWF band struc- 
ture (solid lines) and the one calculated from a simplified 
TB model (filled circles) which includes only the inter-site 
hoppings for which the largest matrix element is larger than 
10 meV (see main text). 



only a very limited number of hopping parameters h'^^ 
between closest neighbors. We therefore identify a min- 
imal subset of hopping parameters, corresponding to in- 
tersite distances |AR|/ac S {1, v^, 2, 3}, i.e. where only 
hopping between sites, for which the leading term (i.e. 
the corresponding matrix element with largest magni- 
tude) is larger than 10 meV, are considered, while the 
rest is set to zero. This model yields an overall very 
good agreement with the full MLWF band structure (see 
inset in Fig.|4]), deviating not more than 0.11 eV for any 
k-point on the 10x10x10 k-point grid used. On the 
other hand, a TB model where only the hopping am- 
plitudes between nearest and next-nearest neighbors are 
taken into account leads to deviations of up to 0.29 eV for 
some k-points, which might still be acceptable for certain 
purposes. However, the overall bandwidth for the latter 
model is reduced by about 0.2 eV compared to the full 
MLWF band structure. 

The calculated matrix elements of the real space ma- 
trix elements for nearest and next nearest neigh- 
bor hopping as well as the corresponding on-site terms 
(AR = 0) are summarized in Table |ll] Here and in 
the following we use the abbreviated notation h^, cor- 
responding to AR = zbocZ, and h^^, corresponding to 
AR = ac(±x±z) (and analogously for all other cartesian 
directions). We note that in the A-AFM case the trans- 
lational equivalence between the two Mn sites within the 
unit cell is broken, and AR = zLucZ is not a lattice vec- 
tor in this case. Nevertheless, in order to simplify the 
notation, we stick to the site-based index and note that 
for A-AFM order a translation along z is equivalent to 
reversing the two spin projections. In the following we 



7 



TABLE II: Calculated values of the on-site, nearest, and next- 
nearest neighbor matrix elements h^^^ (in meV) for FM and 
A-AFM order within structure (i) for the two different spin 
projections. As described in the text, in the A-AFM case all 
matrix elements refer to the Mn site closest to the origin. 



(a) 


FM(t) 


FM(i) 


A-AFM(t) 


A-AFM(i) 


hi. 


14485.7 


17483.9 


14638.9 


17379.3 


h%2 


14484.3 


17483.7 


14541.7 


17443.3 


/if, 
11 


-648.2 


-512.5 




595.0 


22 


9.1 


-9.2 




-8.5 


'hi 


-155.2 


-135.0 


-172.8 


-130.7 


LX 

"■12 


284.5 


217.9 


281.4 


214.5 


ux 
"22 


-483.9 


-386.7 


-488.8 


-389.8 


uxz 
"11 


37.8 


32.0 




34.7 


uxz 
"12 


-34.1 


-26.6 


-29.3 


-29.0 


uxz 
"22 


-1.8 


1.2 




0.5 


"11 


-21.9 


-14.3 


-16.4 


-14.2 


u^y 

"22 


57.9 


47.4 


51.4 


48.3 



always report hopping amplitudes corresponding to hop- 
ping from and to the Mn site at the origin, the corre- 
sponding parameters for all other sites within the unit 
cell follow from symmetry. Similarly, we do not add a 
spin index to the MLWF matrix elements but instead 
discuss each case separately. 

It can be seen that the hopping parameter between two 
|3z^ — r^)-hke MLWFs along the z direction, /if ^ (= t in 
the effective model description), is the leading term for 
the nearest neighbor hopping, and that overall the next 
nearest neighbor hopping is about an order of magni- 
tude smaller than the nearest neighbor hopping. The 
hopping amplitude between two \x'^ — y^)-like functions 
along the z direction, /i^j (= t' in the model description), 
is indeed very small compared to hfi- In the FM case, 
all nearest neighbor hopping amplitudes for the minority 
spin orbitals (except /122) are reduced (to about 75-85%) 
compared to the majority spin channel. This reflects the 
weaker hybridization between minority spin Cg and 0(2p) 
states, leading to more localized minority spin MLWFs 
with reduced hopping amplitudes. For A-AFM order, 
hfi corresponds to the hopping between a local majority 
and a local minority spin orbital, and its value, (92% 
of hfi for FM (t)), is intermediate between the corre- 
sponding FM majority and minority values. The A-AFM 
hopping amplitudes within ferromagnetically ordered x- 
y planes for local majority /minority spin directions are 
very similar to the corresponding FM hoppings (differ- 
ing by less than 5 meV) , with the exception of the (local) 
majority spin hf, value, which is larger than that. Simi- 
lar relations between the FM majority and minority spin 
and A-AFM values are also observed for the next-nearest 
neighbor hoppings. 

It can easily be verified, that the hopping parameters 



for FM majority and minority spin fulfill the relations 
described in Eqs. (|3|-([5|, as required for cubic symme- 
try. However, if the terms proportional to t' = /122 are 
neglected, the corresponding equations are not exactly 
fulfilled. Thus, simply neglecting /122 while keeping all 
other nearest neighbor hopping amplitudes unchanged, 
leads to slight deviations from cubic symmetry. Further- 
more, Eqs. are clearly not fulfilled for the A-AFM 
hopping amplitudes, which reflects the overall tetragonal 
symmetry resulting from the magnetic order. 

This symmetry reduction for the A-AFM case is also 
visible in the on-site matrix elements hi, and /122, which 
differ by about 100 meV. On the other hand, the small 
asymmetry (~ 1 meV) in these on-site terms for FM or- 
der results from small numerical accuracies during the to- 
tal spread minimization (which uses the full k-point grid, 
so that cubic symmetry is not automatically enforced). 

Within the effective two-band model for manganites 
described in Sec. |II A| the Hund's rule coupling leads to 
an on-site spin splitting equal to 2 J (treating Sr as clas- 
sical unit vector). From the calculated on-site MLWF 
matrix elements, we thus obtain a value of J = 1.499 eV 
for the Hund's rule coupling parameter in the FM case, 
and J = 1.370/1.451 eV from the A-AFM on-site terms. 
The differences between these values indicate the limits 
of the assumption of a fixed t2g core spin. We note that 
all these values are slightly larger than the results ob- 
tained in previous LSDA calculations (J = 1.34 eV),'^^ 
which reflects the fact that GGA in g enera l leads to a 
stronger magnetic splitting than LSDA.I^^^^H 

Overall, the results obtained via MLWFs are in a very 
good qualitative agreement with the previous study us- 
ing TB fits to DFT band structures. However, the 
direct comparison between the values calculated from 
MLWF in this work and the values reported in Ref. Hyis 
slightly hampered by the different exchange correlation 
functionals and pseudopotentials used in the two studies. 
The same fitting method as described in Ref . [10] applied 
to the GGA band structure calculated in the present 
work, leads to a nearest neighbor hopping parameter 
t — —688 meV, i.e. slightly larger than the —648 meV 
obtained from the MLWFs. This is due to the larger ma- 
jority spin Cg bandwidth obtained here, — 4.126 eV, 
compared to the value of 3.928 eV reported in Ref. 
Thus, the difference in bandwidth compensates the ne- 
glect of further neighbor hopping in the simple TB fit, 
leading to the apparent very good agreement between 
hfi = —648 meV listed in Table |lT] and the correspond- 
ing value {t — —655 meV) given in Ref. "10^, 

In the following sections, we will analyze the influ- 
ence of the structural distortions only for the on-site and 
nearest-neighbor hopping terms. We have verified that 
the resulting changes in the further neighbor hopping am- 
plitudes do not lead to significant differences in the dis- 
persion characteristics of the eg bands, even though the 
corresponding relative changes of the next-nearest neigh- 
bor hoppings are comparable with those of the nearest- 
neighbor hoppings. 



8 




z r M Az r M Azr 



M A 



FIG. 5: (Color online) DFT band structure (thin lines) for 
the JT-distorted structure (ii): a) majority spin FM, b) mi- 
nority spin FM, and c) A-AFM. MLWF bands are depicted 
as thick/red lines. The Fermi level is indicated by the dashed 
line. 



B. Jahn- Teller distortion — structure (ii) 



As described in Sec. |II C| the staggered JT distortion, 
= ±Qq, leads to a unit cell doubling within the x-y 
plane. In the case of FM order, we therefore construct 
two pairs of Cg MLWFs for each spin channel, localized 
at the two different Mn sites within the unit cell, while 
for A-AFM order we construct four pairs of MLWFs, lo- 
calized at the four different Mn sites within the corre- 
sponding unit cell (for global spin up projection only). 
The same approach for choosing the energy window for 
the disentanglement procedure was used as described in 
the previous section. 

The calculated DFT band structure and e^-hke MLWF 
dispersion for the JT distorted structure (ii) are shown 
in Fig. |5] As a result of the unit cell doubling, there 
are now 4 and 8 bands with Cg character per spin chan- 
nel for the FM and A-AFM order, respectively. As for 
the cubic perovskite structure, the calculated MLWF dis- 
persion largely follows the DFT band structure, except 
where there is strong hybridization with states of a differ- 
ent orbital character. It can be seen that several degen- 
eracies and potential band crossings, which would result 
from a simple "backfolding" of the cubic band-structure 
onto the smaller tetragonal BZ, are lifted due to the JT 
distortion. This can be seen for example for the FM 
majority spin bands, where the highest-lying band along 
FZ acquires some dispersion, leading to a splitting of the 
higher energy eg states at Z. Similarly, the degeneracy of 
the two lowest-lying Cg states at F is lifted, and a poten- 
tial crossing of eg bands is prevented between F and M. 
The latter splitting, together with the reduced dispersion 
along FZ for A-AFM order, appears crucial for the open- 
ing of an energy gap in the JT-distorted A-AFM ordered 
structure (Fig. [sj;). 

To further analyze the influence of the JT distortion 
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FIG. 6: (Color online) MLWF Hamiltonian matrix elements 
ftf'j^ as function of the JT distortion. Large/black and 
small/red symbols correspond to FM and A-AFM order, re- 
spectively. Matrix elements associated with pure (local) ma- 
jority and minority spin character are shown as triangles 
pointing up and down, respectively. Closed circles in (b) rep- 
resent the A-AFM h^-^ hopping. 



on the Cg electronic structure, we perform a series of cal- 
culations where we gradually change the oxygen posi- 
tions from the ideal perovskite structure (i) to the fully 
JT distorted structure (ii), according to Eq. (16 1, and 



monitor the resulting changes in the MLWF Hamilto- 
nian matrix elements. In all these calculations, we use 
the same energy windows of [12.0, 17.5] eV, [15.9, 20.0] eV 
and [12.0, 20.0] eV for the disentanglement in the case 
of FM majority, FM minority and A-AFM, respectively. 
The resulting MLWF matrix elements are depicted in 
Fig. |6] As discussed in the previous section, we report 
only hopping from and to the Mn site at the origin. The 
hopping amplitudes corresponding to other sites in the 
unit cell follow from symmetry. We find a strong lin- 
ear dependence on the JT distortion for both the off- 
diagonal on-site matrix elements (Fig. |6^) as well as 
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for the ofF-diagonal in-plane hopping h^2/2i (^ig-j^/*^)- 
All other on-site and nearest neighbor hopping matrix 
elements show only a weak or moderate quadratic de- 
pendence on Qfii. 

Within the model described in Sec. Ill Al the sole effect 
of the JT distortion (Q^, Qfj^) is a linear coupling to the 
on-site terms at site R according to: 



(17) 



In our case = and = iajiQg; eg is the on-site 
energy of the Cg orbitals. It can be seen from Fig.[6i that 
the off-diagonal element hf{2 indeed shows a linear aepen- 
dence on a, consistent with Eq. (17 1. The correspond- 
ing slope, — AQo = 482 meV, is nearly identical for the 
FM majority and A-AFM local majority spin elements, 
whereas it is significantly smaller for the (local) minority 
spin matrix elements (— AQq = 246/155 meV). This in- 
dicates that the JT splitting is also a ligand- field effect, 
i.e. it is mediated by hybridization with the surround- 
ing oxygen orbitals, which, as pointed out previously, is 
stronger for the energetically lower majority spin states. 
The values for the JT coupling constant A obtained from 
the data shown in Fig. 6a are 3.19 eV/A, 1.63 eV/A, and 
1.02 eV/A, for majority, FM minority, and A-AFM local 
minority spin states, respectively. We note that the value 
of A obtained for majority spin is approximately a fac- 
tor of two larger than the value obtained from the fitting 
procedure described in Ref. [10] As we will discuss in 
more detail below, the source for this discrepancy is the 
strong linear splitting observed for the off-diagonal in- 
plane nearest neighbor hoppings ^-22/21' which is induced 
by the JT distortion (see Fig. |6|:/d). 

This splitting between /i^2/2i *g^in results from the 
underlying hopping between atomic Mn(eg) and 0(p) 
states, which (in leading order) depends linearly on the 
Mn-0 distance. Since this dependence will be different 
for the |3z^ — r^) and \x'^ — y^) orbitals, it can easily 
be verified that the effective hopping across a combina- 
tion of one long and one short Mn-0 bond within the 
X'U plane between two different eg orbitals will also de- 
pend linearly on the JT distortion, whereas the effective 
hopping between the same type of Cg orbitals will show 
only a quadratic dependence. We have verified, by con- 
structing atomic-like Wannier functions for both Mn(eg) 
and 0(p) orbitals (corresponding to larger energy win- 
dows) , that indeed the dependence on the Mn-0 distance 
is much stronger for the hopping amplitude between the 
|3z^ — r^)-type orbital and a neighboring 0{p) orbital 



than for the corresponding 



y )-type hopping, con- 



sistent with the observed splitting in the effective hopping 
ampHtudes ^12/21 shown in Fig.j6j:/d. 

It can be verified within a TB model where the lin- 
ear splitting between /ifj ™d ft-fj^ (and analogously for 
the hopping along the y direction) is taken into account 
via one extra parameter derived from the MLWF data, 
that this splitting partially cancels the effect of the on- 



site JT term on the band dispersion. In particular, the 
JT-induced "gap" between the second and third eg band 
between F and M is reduced by increasing the /ii2/^2i 
splitting, whereas it is enhanced by increasing the JT 
coupling strength A. Thus, the band dispersion result- 
ing from reduced A and no splitting between /i^2 ^^d K^i 
looks very similar to the one obtained from the MLWF 
parameters (i.e. including the spitting between h^2/2'\)- 
This is the reason why the fitting of the DFT band struc- 
ture on a TB model that does not incorporate a /if2/^2i 
splitting (see Ref. 10) leads to a smaller value of A than 
the one obtained from the MLWF parameters. An inter- 
esting question arising from this is whether, despite the 
very similar band dispersion, the two different TB pa- 
rameterizations would lead to noticeable differences in 
calculated ordering temperatures for the collective JT 
distortion. 

The differences between the off-diagonal in-plane hop- 
ping parameters induced by the JT distortion indicate 
changes of the MLWFs themselves, i.e. the JT distortion 
alters the basis-set of a MLWF-based TB model. We note 
that this is an unavoidable result of the effective "two- 
band" picture. The definition of a distortion-independent 
basis-set is only possible within a full d-p TB model, 
based on truly atomic-like functions. On the other hand, 
a splitting between ^nd can in principle also re- 



sult from a unitary mixing of the |3z 



and 



-y 



basis functions. In order to check whether (at least part 
of) the observed splitting is due to such a mixing, we have 
applied a local unitary transformation between the two 
MLWFs on each site, and studied the resulting changes 
in the MLWF matrix elements. In essence, we find that 
it is impossible to retrieve the "cubic symmetry" , i.e. the 
form described in Eqs. (|3|-(|5| and (17 1, simultaneously 
for h*', h^, and h^, and that a transformation of one of 
these terms to the desired form in general increases the 
corresponding deviations in the other two terms. It ap- 
pears that the basis functions resulting directly from the 
maximum localization procedure using initial projections 
on atomic |3z'^ — r^) and \x'^—y'^) functions represent the 
best overall compromise. 



The leading hopping term in z direction, h\-^ (Fig.|6]3), 
exhibits only a weak quadratic change as a function of 
aii. We also find a similar weak quadratic dependence 
on the JT distortion in the hopping parameters h^if22 
(not shown), and a moderately strong quadratic change 
in the on-site diagonal matrix elements (Fig.[6^/f), which 
introduces a splitting of about 150 meV between h^^ and 
/122 for the fully JT distorted structure. 

Finally, we note that the Hund's rule coupling pa- 
rameters derived from the local spin splitting between 
MLWFs obtained for the fully JT distorted structure 
(J = 1.499/1.484 eV for FM order, J = 1.345/1.465 eV 
for A-AFM order) are not significantly changed compared 
with the ones obtained for structure (i). 
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FIG. 7: (Color online) DFT band structure (thin lines) for the 
purely GFO-distorted structure (iii): a) FM majority spin, b) 
FM minority spin, c) A-AFM. MLWF bands are depicted as 
thick/red lines. The Fermi level is indicated by dashed lines. 



C. GdFeOa-type distortion — structure (iii) 

The band dispersion calculated for the purely GFO- 
distorted structure (iii) is presented in Fig. [7] The ro- 
tation and tilting of oxygen octahedra in structure (iii) 
distorts the ideal 180° Mn-O-Mn bond angle, which is ex- 
pected to reduce the hopping amplitudes. Indeed, it can 
be seen in Fig. [7] that the GFO distortion leads to sig- 
nificantly smaller bandwidth (2.951 eV and 2.139 eV for 
FM majority and minority spin, respectively, compared 
to 4.126 eV and 3.156 eV in the undistorted structure 
(i)). As a result, the FM majority spin bands become 
completely separated from the lower-lying t2g bands and 
the La((i) bands at higher energy. Unlike in the JT dis- 
torted structure (ii), the system stays metallic for both 
FM and A-AFM order. 

Since the unit cell for structure (iii) is quadrupled with 
respect to the cubic perovskite structure, there are now 
8 bands with dominant eg character for each spin direc- 
tion. However, due to the tilt/rotation of the oxygen oc- 
tahedra, "Cg-like" orbitals at a certain site can hybridize 
with "t2g-like" orbitals at a neighboring site, leading to 
bands with mixed eg/t2g character.!^ In the FM case 
this does not represent a problem for the disentanglement 
procedure, since the bands with predominant Eg charac- 
ter are separated from the predominantly t2g bands for 
both spin direction. For FM order, we can therefore con- 
struct four pairs of MLWFs, localized at the four different 
sites within the unit cell, by defining appropriate energy 
windows separately for each spin direction. This is not 
possible in the A-AFM case, where the local minority t2g 
bands overlap strongly with the local majority bands 
in the energy region between 14 eV and 16 eV. In this 
case, the standard disentanglement procedure employed 
for structures (i) and (ii), i.e. defining an energy window 
[12.0,21.0] eV and initializing 8 Wannier functions from 
projections on atomic eg orbitals at the various sites, re- 
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FIG. 8: (Color online) Hamiltonian matrix elements in 
the basis of MLWFs as a function of the GFO distortion. 
Large/black and small/red symbols correspond to FM and 
A-AFM order, respectively. Elements associated with purely 
(local) majority and minority spin character are represented 
by triangles pointing up and down, respectively. 



suits in MLWFs with mixed t2g/eg orbital character. In 
particular, the resulting local minority spin MLWFs ex- 
hibit a rather strong t2g character. 

One possible way to overcome this problem would be 
to construct all 20 d-like MLWFs (5 per Mn site), i.e. 
both eg and t2g orbitals. However, the resulting ML- 
WFs still contain some amount of eg/t2g mixing, and 
the corresponding MLWF matrix elements exhibit sys- 
tematic deviations from the results obtained in the pre- 
vious sections, which are derived from a smaller set of 
MLWFs. In the following, we therefore adopt a differ- 
ent strategy to obtain model parameters for the A-AFM 
case, and construct the 4 local majority and 4 local mi- 
nority spin eg-like MLWFs separately, using two different 
energy windows. From this, we obtain the on-site matrix 
elements h*^ as well as the hopping parameters within 
the x-y plane (and of course all further neighbor hopping 
amplitudes within this plane). On the other hand we do 
not obtain the hopping amplitudes between adjacent 
planes in the z direction, which would connect the two 
separate sets of MLWFs. Similar to the purely JT dis- 
torted case, we analyze the effect of the GFO distortion 
on the eg bands by performing calculations with varying 
degree of distortion, i.e. by changing the oxygen positions 
according to Eq. (16). In this case we always adjust the 



energy window for the construction of the MLWFs to the 
actual eg bandwidth corresponding to a particular am. 
We find that the main effect of the GFO distortion 
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FIG. 9: (Color online) DFT band structure (thin lines) for 
structure (iv): a) majority spin FM, b) minority spin FM, 
and c) A-AFM. MLWF are depicted as thick/red lines. Fermi 
level is indicated by dashed line. 



is indeed a systematic reduction of all hopping ampli- 
tudes by ~ 20 — 30%, consistent with what was re- 
ported in Ref. [Tni Fig. [8^ shows the overall reduc- 
tion for all obtained nearest neighbor hopping ampli- 
tudes for both FM and A-AFM order, while Fig. (Sja re- 
solves the reduction factors of the various hopping am- 
plitudes for full GFO distortion {ani = 1). It can be 
seen, that even though there is a significant spread in 
the reduction factors for the various hopping parameters, 
the overall reduction can approximately be described as 
h^/^(Q;iii) = h^/^(0) (l — "Hoi^ii), with an average value of 
7? = 0.26. 

In addition to the changes in the nearest neighbor hop- 
ping amplitudes, we also observe a quadratic decrease 
of the on-site diagonal matrix elements as a function of 
the GFO distortion (Fig. [sji/d), with a similar magni- 
tude for both orbitals and different magnetic order. This 
can be understood again from the underlying hopping 
between atomic p and d orbitals. Since the effective 
eg bands correspond to the antibonding combination of 
these atomic orbitals, a reduction of the underlying p-d 
hopping amplitudes results in a decrease of the F-point 
energy of the Cg states. The Hund's rule coupling param- 
eter J = 1.502 eV obtained from the on-site splitting for 
FM order and = 1 is very similar to the correspond- 
ing value for the cubic perovskite structure. 



D. Combined Jahn- Teller and GdFeOs-type 
distortion — structure (iv) 

So far we have analyzed the individual effects of the 
JT and GFO distortion. We now discuss whether the su- 
perposition of both distortions gives rise to any changes 
in the MLWF matrix elements that go beyond a simple 
superposition of the individual effects. The correspond- 
ing band structure and MLWF dispersion for structure 




FIG. 10: (Color online) MLWF Hamiltonian matrix elements 
as function of combined JT and GFO distortion. Large/black 
and small/red symbols correspond to the FM and A-AFM 
order, respectively. Elements associated with purely (local) 
majority and minority spin character are represented by tri- 
angles pointing up and down, respectively. 



(iv), i.e. the combined JT and GFO distortion, is pre- 
sented in Fig. [9] It can be seen that the band structure 
in this case closely resembles the one of the purely GFO 
distorted structure (iii) , Fig. [t] but with the additional 
JT-induced effects (avoided ba nd-cro ssings and lifted de- 
generacies) as described in Sec. |IIIB Note that, as in the 
purely JT distorted structure, the FM case is metallic, 
whereas a band gap opens only for A-AFM order. 

As described in the previous section we construct 8 
MLWFs per spin direction for the FM case and two sep- 



arate sets of 4 local majority and 4 local minority Cg-like 



MLWFs for the A-AFM case. Fig. [10] shows the evolu- 
tion of selected MLWF matrix elements as a function of 
distortion. The atomic positions are changed according 
to Eq. ( [l6| with x = iv. By comparing Fig. 10 1 with 
Fig. [6^, it can be seen that the GFO distortion does also 
significantly reduce the on-site matrix elements (to 
sa 75 — 80 %) , which are otherwise proportional to the JT 
distortion. This is further evidence for the ligand-field 
nature of the JT coupling, i.e. that it is mediated by 
the Mn-0 hybridization (which is reduced by the GFO 
distortion). Furthermore, it can be seen that the leading 
hopping along z, h\i, follows very closely the trend ob- 
served for the purely GFO distorted structures. In the 
case of the off-diagonal hopping amplitudes within the 
x-y plane, the superposition of GFO-distortion-induced 
reduction and JT-induced splitting leads to an initial in- 
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crease of for small distortion, followed by a decrease 
for larger aiv Overall, the observed trends can indeed 
be well understood as independent superposition of the 
individual effects of JT and GFO distortions. We note 
that the kinks observed in some of the minority spin hop- 
ping terms around ajv ~ 0.4 result from the opening of 
the gap between Cg-like and t2g-likG minority spin bands 
for this amount of distortion, which represents a certain 
"discontinuity" for the disentanglement procedure. 



tional contribution to the in-plane hopping: 



E. Simplified TB models for LaMnOa in the full 
experimental Pbnm structure (v) 



The analysis presented so far showed that the effect of 
different structural distortions on the eg bands can, to 
a good extent, be treated independently of each other. 
In this section, we attempt to incorporate the most sig- 
nificant effects described in the previous sections into a 
refined effective TB model. Then, in order to test the 
accuracy of the resulting parameterization, we compare 
the resulting band dispersion with the full GGA and 
MLWF band-structure, calculated for the full experimen- 
tal Pbnm structure of LaMnOs and A-AFM order. 

For the refined TB model we introduce different hop- 
ping amplitudes for local majority /minority spin pro- 
jections to describe the hopping between ferromagnet- 
ically aligned nearest neighbors within the x-y planes 
(t^^/t^^), and an intermediate value for the nearest neigh- 
bor hopping between antiferromagnetically aligned near- 
est neighbors along the z direction (i^^), i.e. hopping 
between two different local spin projections. This is in 



accordance with our results presented in Sec. Ill A For 
the corresponding hopping amplitudes we use the values 
of /if 2 calculated for the ideal cubic perovskite structure 
(see Table |ll]) for FM and A-AFM order, which are then 
reduced by the same factor (1— Tya?;;), where am describes 
the amount of pure GFO distortion. Apart from these 
modifications we assume the usual cubic symmetry of the 
nearest neighbor hopping matrices, i.e.: 



t'' (±acz) = (l-ryafii)t 



2,,ss' (I 



,0 0/ ' 



(18) 



t-'(±aei) = (1 - ^a?u)t^^' (j^ f] , (19) 



(and analogously for t*'* (±acy)). Note that s and s' in 
these equations should be read as a local spin index, i.e. it 
designates the spin projection relative to the orientation 
of the local core spin. We use the average value ij — 0.26 
determined in Sec. IIII C] 

The JT- induced splitting of the non-diagonal elements 
of the hopping matrix within the x-y plane discussed in 



At(±acx) = AQoQ;ii(l - riani) ( _i 



1 



(20) 



(and analogously for At(±acy)). Here, an describes 
the amplitude of the staggered JT distortion, i.e. 

— iQgttii, and the parameter A is determined from 
the average splitting over all hopping amplitudes in the 
purely JT-distorted structure (shown in Fig. [6j:/d). In 
addition, wc include the usual on-site JT effect in essen- 
tially the same form as described in Eq. but with 
a spin-dependent JT coupling constant that is also re- 
duced by the GFO distortion (with the same factor as 
the hopping amplitudes): 



A ^ A^(l - lya^iO 



(21) 



We note that the orthorhombic strain in the experimental 
structure of LaMnOa gives rise to a homogeneous 
component to the JT distortion, i.e. the same 
on all sites, which we take into account within the model 
by using the same coupling constant A* as for the 
component. 

We also include hopping between next nearest neigh- 
bors and between second nearest neighbors along the 
cartesian coordinate axes in the refined TB model, but 
we do not consider any spin-dependence of the corre- 
sponding hopping amplitudes. We describe the hopping 
between next-nearest neighbors by spin-independent pa- 
rameters t^** corresponding to the hopping between two 



-type orbitals along the ±acX ± UcY directions. 



The parameter t^^ is taken as spin average over the corre- 
sponding MLWF matrix elements h^i calculated for the 
cubic structure. All other hopping matrix elements be- 
tween next nearest neighbors are determined from this 
via the following relations, which are derived assuming 
cubic symmetry and indirect hopping only (see Ref . llOp : 



= t^y (1 - (J f 

2 N / 1 



(1 - T]a 



-3 



(22) 
(23) 



Sec. IIIB is incorporated in the TB model as an addi- 



The same GFO-distortion-induced reduction as for the 
nearest neighbor hopping matrices is applied. The hop- 
ping between second nearest neighbors along the coordi- 
nate axes [t(±2acx), t(±2acy), t(±2acz)] is included ac- 
cording to the ideal cubic symmetry relations described 
by Eqs. ([3|-([5|, with Cc replaced by 2ac, t' = 0, and 
t — t^^, where is estimated from the MLWF matrix 
elements for the purely GFO distorted structure. We 
note that the reduction of this parameter compared to 
the undistorted case is significantly stronger than for the 
nearest (and next nearest) neighbor hopping amplitudes. 
Furthermore, the hopping between third nearest neigh- 
bors along the coordinate axes [t(±3acx), t(±3acy), 
t(±3acZ;)], that was considered in Sec. Ill A| becomes 




negligible as result of the GFO distortion. 

Finally, we include the Hund's rule coupling in the re- 
fined TB model using the standard form [Eq. Q] with 
an average value of J obtained from the MLWF on-site 
splitting. In order to relate the obtained TB bands to 
the full GGA and MLWF band-structures, we determine 
the on-site energy eo as the spin and orbital average of 
the corresponding matrix elements for the A-AFM 
experimental structure. 

The values of all parameters used in the refined TB 
model are summarized in Table Hill The JT distortion in 
the experimental Pbnm structure corresponds to = 
1, Qg = -0.161 A, and = -0.048 A, and the corre- 
sponding amphtude of the GFO distortion is am — 1 . 

We also compare with a very simple TB model that 
includes only nearest neighbor hopping according to 
Eqs. with t' = 0, and the standard JT and Hund's 

rule coupling as described by Eqs. Q and (|6|. The pa- 
rameters for this model are chosen via typical simpli- 
fied fitting procedures: the nearest neighbor hopping pa- 
rameter —t is obtained as one sixth of the majority spin 
bandwidth W for the fully GFO distorted structure (iii) 
and FM order; the JT coupling constant A is taken from 
Ref. TIT, where it was obtained by fitting a similar TB 
model (including also next nearest neighbor hopping) to 
a DFT band-structure; J is calculated from the spin split- 
ting between FM majority and minority bands at the F- 
point for the cubic structure (i); and cq is fitted such 
that the Fermi energy is aligned with the DFT calcula- 
tion value. 

Figure |11^ shows the band dispersion obtained from 
the GGA calculation for the full experimental struc- 
ture (v) and A-AFM order as well as the correspond- 
ing MLWF bands. Fig. [TTja/c shows the comparison be- 
tween the MLWF bands and the two different simplified 
TB models. It can be seen that the orthorhombic lattice 
strain and La displacements do not lead to significant 
qualitative changes in the band-structure as compared to 
structure (iv) (see Fig.|9}:). The comparison between the 
MLWF dispersion and the refined TB model (Fig. [llj)) 
shows that, despite the many simplifications made, this 



FIG. 11: (Color online) (a) DFT bands (thin/grey lines) and 
MLWF bands (thick/red lines) for the A-AFM experimental 
Pbnm structure (v). Comparison of the MLWF bands with 
the refined TB model (b) and with the simple nearest neighbor 
TB model (c) The Fermi level is indicated by dashed lines. 

model reproduces the MLWF bands to a remarkable ac- 
curacy. The only major discrepancy can be seen for the 
lowest-lying local minority band along F-Z at ~ 16 eV, 
which is slightly lower than the corresponding MLWF 
band. This can be traced back to an overestimation of 
the /if2 hopping amplitude, which results from the fact 
that we use the same reduction factor ry for all hoppings. 
As can be seen in Fig. [SJd, /ifj is affected more strongly 
by the GFO distortion than any other nearest neighbor 
hopping (for A-AFM order). The very simple nearest 
neighbor TB model depicted in Fig. [TT] ; deviates much 
stronger from the MLWF band structure than the refined 
model, but still captures the overall dispersion surpris- 
ingly well. Consistent with our analysis from the previ- 
ous sections, the deviations are more pronounced for the 
energetically higher local minority spin bands, which is 
clearly due to the neglected spin dependence of the hop- 
ping. As discussed in Sec. |IIIB| the smaller JT couphng 
constant used in the simple model partially cancels the 
missing effect of the JT distortion on the inter-orbital in- 
plane hopping parameters, leading to the relative good 
agreement of the simpler model with the MLWF bands 
around the Fermi level. 



IV. SUMMARY AND CONCLUSIONS 

We have shown that the construction of maximally lo- 
calized Wannier functions together with the disentangle- 
ment procedure described in Ref. T6 can be used to ex- 
tract effective eg bands in LaMnOs even for cases where 
these bands are strongly entangled with other states. 
This procedure thus provides a very robust way for ex- 
tracting the "correlated subspace" used for example in 
DFT-I-DMFT calculations. 

We have used this procedure to obtain a TB param- 
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eterization of the Cg bands for different structural mod- 
ifications of LaMnOs with both FM and A-AFM order. 
By monitoring the effect of the individual distortions on 
the MLWF matrix elements, we can assess the quality 
of the various approximations and simplifications that 
are commonly used in model Hamiltonians for mangan- 
ite systems. In particular, we find the following: 

• While the nearest neighbor hopping is clearly domi- 
nant, the further neighbor hopping along the carte- 
sian axes decays rather slowly. Taking into account 
nearest, next-nearest, as well as second and third 
nearest hopping along the cartesian axes leads to 
deviations of less than 0.11 eV from the (cubic FM) 
DFT band structure. 

• In addition to the linear on-site coupling to the 
JT distortion, we observe a strong effect on the 
in-plane hopping amplitudes between different or- 
bitals. The corresponding splitting, which is due to 
the underlying Mn-0 hopping, partially cancels the 
effect of the on-site term on the band dispersion, 
which has a strong influence on the determination 
of the local JT coupling strength. 

• The GFO distortion leads to an overall reduction 
of all hopping amplitudes by about 25-30 %, and 
also reduces the local JT splitting. This reduction 
is due to the weaker hybridization between Mn(eg) 
and 0(p) states for non-180° bond angle. 

• The higher energy of the (local) minority spin 
states reduces the hybridization between the cor- 
responding atomic eg and 0{p) states, leading to 
reduced hopping amplitudes and JT coupling com- 
pared to the majority spin states. 

• The splitting between (local) majority and minority 
spin states is generally well described by the local 
Hund's rule coupling, even though small variations 
in the corresponding J values indicate the limits of 
the core spin approximation. 

It is apparent that the most crucial deviations from the 
simple two band description are a result of the underly- 
ing Mn-0 hybridization. Nevertheless, we have shown 



that a refined TB model that incorporates the effects de- 
scribed above using the parameters listed in Table III 
reproduces the DFT band structure calculated for the 
full experimental crystal structure of LaMnOa with re- 
markable accuracy. Whether this accuracy, at the prize 
of more parameters in the model, is desirable depends of 
course on the specific application of the model descrip- 
tion. 

Furthermore, our analysis shows that the effects of the 
various distinct structural distortions present in LaMnOa 
are (to a good approximation) independent from each 
other and can therefore be assessed individually. How- 
ever, the GFO distortion has to be taken into account to 
obtain the correct magnitude of the Jahn- Teller coupling. 

In comparison with the manual TB fits presented in 
Ref. ^10, the construction of MLWFs is less biased and 
more universally applicable. It allows to calculate pa- 
rameters of the model instead of fitting them to either 
experimental or computational data. In particular, it is 
possible to obtain accurate TB representations even for 
rather complex band structures. However, care has to be 
applied when parameters corresponding to a more com- 
plex parameterization are used for simpler models. For 
example, using the MLWF nearest neighbor hopping am- 
plitudes within a simple model that neglects all further 
neighbor hoppings, can lead to a significant underesti- 
mation of the total bandwidth, so that in certain cases 
a parameterization with renormalized nearest neighbor 
hoppings, leading to a more accurate total bandwidth, 
might be more desirable. The analysis presented in this 
work demonstrates that, depending on the specific ap- 
plication at hand, MLWFs can in principle be used to 
construct more and more refined TB parameterizations 
which lead to realistic, materials-specific band structures 
with very high accuracy. 
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